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We compare pseudopure state ensemble implementations, quantified by their initial polarization 
and ensemble size, of Grover's search algorithm to probabilistic classical sequential search algo- 
rithms in terms of their success and failure probabilities. We propose a criterion for quantifying 
the resources used by the ensemble implementation via the aggregate number of oracle invocations 
across the entire ensemble and use this as a basis for comparison with classical search algorithms. 
We determine bounds for a critical polarization such that the ensemble algorithm succeeds with a 
greater probability than the probabilistic classical sequential search. Our results indicate that the 
critical polarization scales as TV -1 / 4 where N is the database size and that for typical room temper- 
ature solution state NMR, the polarization is such that the ensemble implementation of Grover's 
algorithm would be advantageous for N > 10 22 . 

PACS numbers: 03.67.Lx 



I. INTRODUCTION 



Conventional implementations of quantum algorithms 
entail application of unitary transformations and pro- 
jective measurements to a single, multi-qubit quantum 
system which is initially prepared in a known pure 
state [HQ,|3, 4]. In the alternative ensemble paradigm, of 
which solution state nuclear magnetic resonance (NMR) 
is the most prominent example, the algorithm is imple- 
mented on an ensemble of multiple identical, noninter- 
acting quantum systems d S, 0, S, S El El El El E| • 

The principal differences between the two paradigms 
manifest themselves in the preparation, measurement 
and post-measurement stages and arise from: (i) the im- 
possibility of isolating any given ensemble members and 
applying preparation, evolution and measurement oper- 
ations to these alone while ignoring the remaining en- 
semble members and (ii) the fact that available initial 
states are mixed. The principal similarity between the 
paradigms is that it is possible to apply identical uni- 
tary evolution operators to each ensemble member at any 
stage of the algorithm. A variety of schemes have been 
developed to deal with the initialization, measurement 
and algorithm output issues in ensemble imp lement at ions 
of quantum algorithms @, 0, El EI El E| ■ Collectively 
these issues are resolved via the use of a pseudopure ini- 
tial state^for who se prep aration there are several known 
schemes 0, @, El 13 E| , described by the density oper- 
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where n is the number of qubits, |^) is the initial pure 
state required by the conventional (single quantum sys- 
tem) implementation of the algorithm and ^ e 1 is 
called the polarization. Under the unitary transformation 
for the algorithm, E/ a ig, the density operator transforms 
to 
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where 



final/ 
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is the final, pre-measurement state of the quantum sys- 
tem in the conventional implementation of the algorithm. 
This gives the appearance of a conventional realization 
of the algorithm on a subset of the ensemble which is ini- 
tially described by the density operator \ipi)(tpi\, which 
corresponds to a pure state. In most proposals for and 
realizations of ensemble implementations the algorithm 
output is extracted from the expectation values of single 
qubit traceless observables @, 0, El El > as we describe 
in Sec. Mil Such expectation values are deterministic. 
However, in real implementations the number of ensem- 
ble members is finite and expectation values can only 
be approximated by sample averages. Thus, regardless 
of the nature of the original quantum algorithm, ensem- 
ble implementations are probabilistic (except possibly in 
the special case where e = 1 which reduces to a conven- 
tional implementation). Although this basic feature has 
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been noted @, i there has been limited discussion of 
its implications with regard to reasonable comparisons 
of ensemble implementations of quantum algorithms to 
probabilistic classical competitors. Such issues have been 
addressed in the context of various ensemble implementa- 
tions of single bit output algorithms such as the Deutsch- 
Jozsa algorithm ppl . |2l| , resulting in bounds of the polar- 
ization beneath which a classical probabilistic algorithm 
will succeed with greater probability than the ensemble 
implementation of the quantum algorithm. 

In this article we compare the performance of an en- 
semble implementation of the Grover search algorithm to 
probabilistic classical sequential searches. Our aim is to 
determine the polarization as a function of both ensemble 
and database size such that the ensemble implementation 
of the quantum algorithm succeeds with greater probabil- 
ity than its classical competitor. The primary distinction 
between this and previous work on the Deutsch-Jozsa al- 
gorithm is that the output of the Grover algorithm is de- 
termined from measurements on many qubits and that 
these measurement outcomes are correlated. 

The article is organized as follows. We describe the 
conventional Grover search algorithm in Sec. [II] In 
Sec. Mil we describe an ensemble implementation of the 
Grover search algorithm and explicitly provide a proto- 
col for extracting an algorithm output. The bulk of our 
work is found in Sec. IIV1 which covers the probability 
with which ensemble implementations succeed and uses 
these as a basis for comparisons with classical probabilis- 
tic searches. Finally the appendices contain elaborations 
of the mathematics behind our results. 



1, . . . , 2" — 1} by setting f(x) = at these points. For 
this reason we shall restrict discussion of the quantum al- 
gorithm to cases where N = 2" for some integer n. The 
algorithm unfolds by initializing the qubits to the state 
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followed by repeated applications of the Grover iterate, 
G := DUf, where the unitary "inversion about the aver- 
age" operation is 
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with (c) = ^2^Zq c x /N. Each application of G includes 
one oracle query. For a database containing only one 
marked item, to which case we restrict our considera- 
tion, a standard analysis [l], [23j] demonstrates that after 
q applications of G the state of the system is 



where 



P„ := cos 
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II. GROVER'S SEARCH ALGORITHM: 
CONVENTIONAL IMPLEMENTATION 

The simplest unstructured database search problem in- 
volves a database that conceals a single marked item in 
one of N possible locations, conveniently represented by 
the integers {0, 1, . . . N — 1}. The task is to find the lo- 
cation of the marked item, s, with the help of an oracle, 



/(*) == 



if x s 

1 if x = s. 



(4) 



A classical sequential search proceeds by evaluating the 
oracle at distinct database locations and on average re- 
quires N/2 such oracle queries to determine the marked 
item's location. 

The Grover search algorithm [l|, [52] requires n = 
|~log 2 N~\ qubits whose computational basis states 
{\x), x = 0, . . . 2™ — 1} can be regarded as representing 
database locations. The oracle is implemented via a uni- 
tary transformation, defined on the computational basis 
states as 



U f \x) := 



(5) 



and extended linearly to all superpositions. The ora- 
cle is easily extended to the database locations {N, N + 



and < 8 < ir depends on TV via 
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The probability that the final computational basis mea- 
surement will yield s is 



Pr (correct) = |(s|^i)| = sin 



2 _ . 2/2^+1 



(12) 



and the measurement returns the marked item's location 
with certainty whenever (2q + 1)0/2 is an odd multi- 
ple of tt/2. The lowest value of q which ensures this is 
it/ (26) — 1/2 but for an arbitrary database size, 9 may 
not be such that this is integral. The standard proto- 
col for the Grover algorithm requires that the number of 
applications of G is 



(13) 



where nint (x) is the nearest integer to x. The prob- 
ability with which a computational basis measurement 
will yield s correctly can be ascertained by noting that 
7r/(20) — 1 ^ q s td ^ 7r /(26'). At either extreme, the prob- 
ability that the computational basis measurement yields 
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s is cos 2 (0/2) = 1 — 1/ N; thus the probability of success- 
fully locating the marked item after q si( x oracle invoca- 
tions is at least 1 — 1/TV. It can be shown that, throughout 
this range, the probability that the final computational 
basis measurement yields s is bounded from below by the 
probability at the extremes of the range. Typically, the 
search problem is considered for TV ^> 1, in which case 
the Grover algorithm yields the marked item's location 
with near certainty. The number of oracle invocations re- 
quired for the standard algorithm can be approximated 
from Eq. (JTTJ) 



III. GROVER SEARCH ALGORITHM ON A 
ENSEMBLE OF QUANTUM SYSTEMS 



and thus 



q st d = runt 
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We consider an ensemble implementation of Grover's 
algorithm starting with a pscudopure state of the type 
given by Eq. (!]). Eqs. ((2|) and JH]) imply that the density 
operator after q applications of G and prior to measure- 
ment is 
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The standard protocol for determining the algorithm outcome, and one which eliminates the term proportional to /, is 
based on the expectation values of traceless single qubit observables [H, [lj|. For the j th qubit, single qubit expectation 
values can be computed from Pcl^red* the final reduced density operator for the qubit. For Grover's algorithm, 
Eq. (fT6|) implies that, after q applications of G, 



,0) 



1 



final red 2 



l-e 



\a a \ 2 N- 1 



TV - 1 



J + e 



K| 2 TV-1 



TV 

1 

2 

Vtv^"T 



i-KI 2 

TV- 1 



TV - 1 

\sj)(s j ®l\+\sj®l){s j \ 



Ctnd* 



(17) 



where Sj is the j bit of s. The expectation value, 

<^> W =Tr(a^ alred ) 

, ... |a | 2 TV-l , , 

= (_!)'» g l (18) 

depends on Sj . Thus Sj can be determined whenever the 
expectation value of a z can be measured for the j th qubit 
and the measurement resolution allows for distinction be- 
tween -K(|a 9 | 2 TV-l)/(TV-l) and-e(|a 9 | 2 TV-l)/(TV- 
1). This opens the possibility for a truncated version of 
Grover's algorithm, which uses fewer oracle invocations 
(i.e. q < q st d) than required by the standard version of 
Grover's algorithm, subject to the proviso that \a q \ is 
large enough to distinguish between the two possible ex- 
pectation values [TH . If these requirements are satisfied, 
then an ensemble implementation based on expectation 



values will be deterministic in the sense that a given value 
of s will always yield the same outcomes and the value 
of s can be determined with certainty. 

For an ensemble containing a finite number of mem- 
bers, M, precise expectation value measurements are ide- 
alizations. Their purpose is to suggest protocols that 
produce algorithm outputs based on appropriate sample 
averages of projective measurement outcomes for individ- 
ual ensemble members. The key notion is that, for a suf- 
ficiently large M, sample averages of projective measure- 
ment outcomes give close approximations to expectation 
values. The typical protocol has already been described 
for single bit output algorithms [2(| and is readily ex- 
tended to multiple bit output algorithms. First, a compu- 
tational basis measurement is performed across all qubits 
for each ensemble member. Consider the jth qubit of en- 
semble member k (where 1 ^ k ^ M). The projectors 
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for the two possible computational basis measurement 



outcomes arc P Sj := \sj)(sj\ and 



© l)(Sj ® 1| 



and the corresponding measurement outcome are scaled 
to z]jp = +(-l) Sj and = -(-1) s j respectively; this 
translates between the computational basis measurement 
outcomes and the eigenvalues of a z . Second, the sample 
average for the jth qubit is computed 
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and this approximates {a z )^\ Finally, if z^ > (or 
conversely z^> < 0) then this suggests that (u z )^ > 0, 
(or (<J Z )^ < 0) which, according to Eq. I|18p implies 
Sj = (or Sj = 1). Thus a candidate, c := c n . . . c\, for 
the marked item 's location is assigned via: 



Cj =0 



Zj = with probability 1/2, 
Zj = 1 with probability 1/2 and 

= 1. 



> 

z« = 

zW> < 

(20) 

The protocol of Eqs. (fT9|) and (f20|) is equivalent to tak- 
ing a bitwise majority vote of the computational basis 
measurement outcomes over the entire ensemble. If, for 
a given bit, more ensemble members return +1 than —1 
after computational basis measurements, then the intu- 
ition is that Sj = since the outcome (— l) Sj is more 
likely than — (— l) Sj . The string c := c„ . . . c\ constructed 
in this fashion constitutes the ensemble algorithm's out- 
put. 

The probabilities with which single qubit outcomes oc- 
cur are readily determined from the reduced density op- 
erator of Eq. (|17[) . Specifically the probability that an 
ensemble member yields the "correct" value is 



1 + £cff 



where the effective polarization is 



£ e ff := e 



\a q \ 2 N- 1 
TV — 1 
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(22) 



Similarly the probability that an ensemble member yields 
the "incorrect" value is 



= ) = 



1 — Ecff 



(23) 



Unless both e = 1 and \a q \ 2 = 1, the protocol of 
Eqs. fT9"]) and (|2T)|) can give Cj ^ Sj, thus causing the 
algorithm to fail to locate the marked item correctly. In 
this sense the algorithm is probabilistic. Even for perfect 
polarization, this is generally unavoidable since, for the 
standard number of invocations of the Grovcr iterate, i.e. 



9 = fed, l^gl 2 7^ I- However, for the standard number of 
iterations, 



giving 



^ £eff ^ 



N- 1 



(24) 



(25) 



and as N — > oo, e e ff ~ £• Thus for the typical scenario, 
for which N 3> 1 and q = q st d in ensemble implementa- 
tions the probabilistic aspects will be determined pre- 
dominantly by the polarization. We shall henceforth 
only consider the typical scenario, taking e e g = e and 
q = nVN/4. 

The protocol of Eq. ([20]) is one in which it is assumed 
that it is possible to distinguish a discrepancy as small as 
1 between the numbers of projective measurement out- 
comes returning +1 and —1 respectively. This is the best 
resolution case and one can generalize the situation to 
consider cases where the resolution is worse as has been 
done for single-bit output algorithms [20l ]. In this arti- 
cle we shall restrict the discussion to the best resolution 
case. 



IV. STATISTICAL PERFORMANCE OF THE 
ENSEMBLE QUANTUM SEARCH ALGORITHM 

The probabilistic nature of quantum algorithms on en- 
sembles motivates performance comparisons, based on 
the probability with which each algorithm fails or suc- 
ceeds to correctly identify the marked item's location, 
with probabilistic classical competitors. Of the possible 
scenarios, we shall consider the probability of correctly 
identifying every bit s\, . . . , s n of the marked item's lo- 
cation. Eqs ([2T|l and |22|) suggest that this will depend 
on the database size N and the polarization e (for the 
typical scenario q is determined by TV). Also, as with 
any sample average, it is to be expected that the fail- 
ure or success probability will depend on the ensemble 
size M. Denote the probability with which the ensem- 
ble quantum algorithm correctly identifies every bit of s 
by -Psuccoss aii q ( e ' N) and the probability with which it 
fails (i.e. at least one bit of s is incorrectly identified) 
by Pfca all q (^ M, N) = 1— Psucc _ all q ( £ , M, N). Methods 
for calculating and approximating these have been pro- 
vided for ensemble quantum algorithms with single bit 
outputs [2(|. The key difference here is that Grover's al- 
gorithm's output involves multiple qubits and that mea- 
surement outcomes for distinct qubits in an individual 
ensemble member are correlated. Thus, although we can 
still use the methods for single bit outputs [20[ to de- 
termine the probability with which any single bit value 
will be correctly identified, the fact that the measure- 
ment outcomes for different qubits are not statistically 
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independent means that success probabilities for multi- 
ple output bit cases will not necessarily be simple com- 
binations of those for the single output bit case. 

To illustrate this point, consider the case where n = 2 
(i.e. N = 4) and the marked item is located at s = 11 
(in binary notation). Eq. (fTT|) implies that = it/ 3 giv- 
ing q std = 1. Then Eq. © implies \%jj f ) = \s) = |11). 
For an ensemble realization, Eq. (f!Z2]) gives £ e ff = £• This 
yields the joint probability distribution for various mea- 
surement outcomes on both qubits of a single ensemble 
member of Table fl] 

The marginal probability distributions for single qubit 
outcomes are identical; for the k th ensemble member and 
the j th bit, Pr {z^ ] = 0) = (1 - e)/2 and Pr (zjf* = 1) = 
(1 + e)/2. Thus the probability distributions for the in- 
dividual qubit outcomes are not independent except in 
the trivial cases of e = or e = 1. Let hi denote the 
number of times that the outcome for qubit 1 is 1 and 
the number of times that it is 1 for qubit 2. When the 
number of ensemble members is odd the ensemble algo- 



rithm succeeds provided that fci > M/2 and k 2 > M/2. 
It is shown in appendix [5] that the probability with which 
single qubit measurements yield fci and k 2 is 



Bit 2 Outcome 


Bit 1 Outcome 


Probability 








1-e 








4 





1 


1-e 








4 







1-e 


1 






4 


1 


1 


1 + 3e 
4 



TABLE I: Joint probability distribution for N = 4 and s — 11 
(binary). 



Pr(fe,fci) = 



4 M 



min (fei ,&2) 

(M -h- k 2 +l,h- I, k 2 - I, l)\ (1 - e) M ~ l (1 + 3e)' 

=max (0,fei+fc2-Af) 



(26) 



where (ni, ri2, ns, 714)! 
714)! / 1 (n\\n 2 \ny,\ni\) . Then 



M 



Success a ii q (£,M,JV)= £ 



("1 



E 



"2 + n 3 



Pr(fea,Ai). (27) 



7. M + l i M + l 

For algorithms with single bit outputs, the expression 
analogous to Eq. ([??]) is a cumulative binomial distri- 
bution, which can be expressed as an incomplete beta 
function, greatly simplifying various calculations and in- 
ferences [20|. We are unaware of comparable techniques 
for re-expressing Eq. (|27[) in cases involving more than 
one output bit. 

Nevertheless certain general conclusions are possible. 
If e = 0, the ensemble is in the maximally mixed state 
throughout the algorithm. Thus for a given qubit, a com- 
putational basis measurement yields either outcome with 
equal probability. The procedure for deciding the algo- 
rithm output amounts to unbiased guessing on each bit 
and these outcomes are no longer correlated. It follows 
that 



1 



1 



Success ali q (0,M,^) = ^ = ^- 



(28) 



IV'finalXV': 



At the other extreme e = 1 and p nna i 
The fact that in the typical scenario \a 
cssarily 1 implies that in general, e c ff 7^ 1, and the mea- 



Qstd I 



final I • 

is not nec- 



surcmcnt outcomes are correlated. However, e e g w 1 and 



success 



aai q (l) M, N) « 1. 



(29) 



for N 3> 1. A more precise lower bound on 

PsucccssanHM-^^O can be obtained via Eqs. (J33D 
and PO]): this implies the result of Eq. (gSJ) for N > 1. 

The absence of techniques for computing or approxi- 
mating P success al i q (e, M, N) for more than a single out- 
put bit motivates the search for bounds on this quantity 
or, equivalcntly and more conveniently, p fail a n q (e, M, N). 
The strategy is to bound this in terms of the single 
bit failure probability, i.e. the marginal probability with 
which the algorithm fails to identify a given single qubit 
correctly regardless of the output for the other bits. As- 
suming, without loss of generality that the marked item 
is located at s = 11... 11, Eq. fl2|) implies that the 
joint probability distribution is symmetrical under in- 
terchanges of qubits. Thus the single bit failure prob- 
ability is identical for all qubits and will be denoted 
p i ono q (e, M, N). The bounding relationship, central to 
the rest of this article, between the single bit failure prob- 
ability and the probability with which the algorithm fails 



Pfail c 



<Pfailali q < " Pfa 



(30) 



where the arguments, (s,M,N), of the probabilities are 
identical and have been omitted for convenience (sec ap- 
pendix [B] for a proof). 
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Following; methods established for single bit output al- 
gorithms [20 ]. 



Pfailonc^^) 
1 M 

\ E 



k=M„ 



1 

+ 2 ^ V 

fe=M-M min +l 



E 



MA ( 1 - e cff \ k f 1 + e cff x 



where £ e ff is given by Eq. (I22[) and 



(31) 



(32) 



The single bit failure probability is a monotonically de- 
creasing function of e [2(| with p fail one q (0, M, iV) = 1/2 
and 



2 2(AT-1) 



ttM 2N - 3 V N - 1 



Af/2 



(33) 

whenever M ^> 1 (see appendix [C]). Thus for A^ ^> 
1, p f ;l one q (l, M, TV) w 0. For fixed polarization e and 
odd M, P ianonc q (£,M,N) = p lailone q (e,M + 1,JV) and 

Pfea onc q (£> M, iV ) r %ftaii onc q (^ Af + 2, AT) with equality 
only if £ off = 0, 1 [20(. 

The exact expression of Eq. (f3~Tj) is not conducive to 
rapid numerical evaluation and is less useful for assessing 
the behavior of p fail ono q (e, M, JV) than certain alterna- 
tives. It can easily be verified by repeated integration by 
parts that 



k—ra 



£ (")/(l-pr- fe = /p(m,n-m+l) 



(34) 



where 



I p (x,y):=^yl f P t»-i(i _ t )«-i (35) 

r(x)r(y) 7 

is the incomplete beta function. Applied to Eq. (|31 [) this 
yields an alternative exact expression, 

Pm onc q (^ Af, JV) =\ [l p (M min , M - M min + 1) 

+ I p (M-M min + l,M min )] 

(36) 

where p = (1 - e)/2. For M > 1, Eqs. ([35]) and ((36]) can 
be combined (see appendix [D]) to give the approximation 



This is closely related to the standard error function and 
is more suitable for numerical computation than the ex- 
act expression of Eq. (|3"Tj) . 



A. Statistical comparison of the ensemble quantum 
and classical sequential search algorithms 

The ensemble quantum algorithm is probabilistic and 
it is reasonable to compare it in terms of success prob- 
abilities to a probabilistic classical sequential search. A 
deterministic classical sequential search proceeds by eval- 
uating the oracle, f(x), at distinct database locations x, 
terminating whenever f(x) = 1. We consider a proba- 
bilistic sequential search in which the oracle is evaluated 
at Q distinct randomly chosen database locations. The 
probability with which this fails to correctly identify the 
marked item is 



Pfai lall C (Q^) = l 



Q 

N 



(38) 



provided that Q ^ N—2 since a search with N—l queries 
will yield the marked item's location with certainty. We 
shall compare the two types of searches by requiring that 
each use the same resources measured in terms of the ag- 
gregate number of oracle invocations (we briefly discuss 
other possibilities later). Thus we regard an ensemble 
implementation invoking q oracle invocations, each on 
M ensemble members, as using equivalent resources to 
a classical sequential with Q = qM oracle queries. The 
central issue in this article is to compare the two types 
of searches provided that Q = qM and to determine con- 
ditions under which the ensemble implementation of the 
quantum algorithm outperforms the probabilistic classi- 
cal sequential search, i.e. such that 



PfaiiaiA^AO < P{ailaU c (qM,N). 



(39) 



The critical polarization e c , which demarcates the bound- 
ary between the two types of algorithm, is defined implic- 
itly via 



P faiia n <l (^M,N)=p iaiiaU c ( q M,N). 



(40) 



In general e c can be expected to depend on M and N 
and will be denoted e c (M, N). Eqs. (28]), j29]) and fiOjl 
imply that a critical polarization in the range e c $J 1 
can only exist provided that 



1 sC qM ^ N. 



(41) 



The left bound is trivially satisfied while the right bound 
implies an upper limit on the ensemble size, 



Mr, 



(42) 



For M > M max , p fail all c (qM, N) = and the classical 
algorithm always outperforms the quantum algorithm, 
except in some (but not all) situations for which e = 
1. For M < M max the critical polarization exists since 
p f ., aU q (e, M, N) is a continuous function of e and the 
classical failure probability lies within the bounds for the 
quantum failure probability. 
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The difficulties in calculating p succoss a i q ( e : M, N) are 
equivalent to those in calculating p { ., al i q (£, M, N) and 
we resort to calculating bounds, originating from the in- 
equality of Eq. (|30p. for the critical polarization. We shall 
show that 

e c nocc (Af, N) < e c (M, N) < e c 8U ff(M, TV) (43) 
where the necessary critical polarization, e c necc , satisfies 



Q = qM imply that 



l (e cnecc ,M > N)=p faiiaU c (qM,N). 



(44) 



and the sufficient critical polarization, £ C suffi satisfies 

Pf ail onc q (£c S u ff ,Af,AO = ^ Pmiall c (qM,N). (45) 

Each of these depend on M and N and they will some- 
times be denoted e cnecc (M,N) and e csu s(M, N). As 
the necessary and sufficient critical polarizations are de- 
fined via the single bit failure probability existing tech- 
niques (2(j| can be used to demonstrate their properties 
and calculate them. 

The existence of the necessary critical polarization fol- 
lows from the fact [2(| that p fail onc q (e, M, N) is a mono- 
tonically decreasing function of e ranging from to 1/2. 
Thus a unique value of the necessary critical polariza- 
tion exists whenever \N/2q] ^ M ^ M max . Note that in 
the typical scenario this implies M max /2 ^ M ^ M max . 
For M ^ \N/2q], the classical sequential search fails 
with probability larger than 1/2 and thus the ensem- 
ble quantum algorithm always outperforms it, giving 
£c necc (A/, N) = 0. Similarly the sufficient critical polar- 
ization exists provided that (1 — n/2)N/q ^ M ^ M max . 
The left hand inequality is satisfied whenever n ^ 2 as 
it does in the typical scenario. Here a unique sufficient 
polarization exists when M Af max . 

The term "necessary" stems from considering e < 
£cnocc, in which case the monotonic behavior of 
Pfaii onc q ( £ > M > N ) witn respect to e and Eqs (J30J) and ((44|) 
imply 



Pfailan q (^^^)>Pfailan C (^^); 



(46) 



here the classical algorithm succeeds with greater prob- 
ability than the quantum ensemble algorithm. In this 
sense, £ c noC c establishes a lower bound on e c . 

The term "sufficient" arises by considering, e > e c su ff , 
for which a similar argument implies 



Pu aa n q (e,M,N)<p {ailall c (qM,N); 



(47) 



here the quantum ensemble algorithm succeeds with 
greater probability than the classical sequential search 
using the same resources. In this regard £ c su ff is an up- 
per bound on e c . This establishes Eq. ([4*3")) . 

The task now becomes one of computing the two 
bounding critical polarizations. Eqs. (j3"8")) and (|4"4"|) with 



Pfail onc q ( £ c necc, M,N) = 1~^- 

ttM 



1 - 



4VN 



(48) 



For fixed N, Eq. (|48|) can be solved numerically for 
£ c nocc(M, N). Data, computed via Eq. ([36]) using MATH- 
EMATICA's implementation of the incomplete beta 
function, is plotted in Fig. Q] 



£(1(T 2 ) ' 



a) 2 



2.5 5.0 7.5 10.0 12.5 M(10 4 ) 



MKT 3 ) 

3 - 

b) 2 



2.5 5.0 7.5 10.0 12.5 M(10 6 ) 

FIG. 1: Critical polarization for a) N = 10 10 and b) N = 10 14 . 
Open squares indicate sufficient polarization and solid dots 
necessary polarization. Note the distinct horizontal and ver- 
tical scales. The rightmost point in each plot corresponds to 
M — M max for the relevant value of N. The critical polariza- 
tion, e c , lies between the two curves. 

In general it can be shown (see appendix |E| that 
£c nccc(Af, N) is a monotonically increasing function of M 
with £ c nGCC — * 1 as M — > M max . A more important con- 
sideration is the scaling behavior of e c necc (M, N) with 
respect to N while fixing M suitably. The similarity in 
the data generated for the two values of N illustrated in 
Fig. Q] suggests an approach based on £ c nGCC as a func- 
tion M/M max for different values of N. In fact, it is easily 
shown that 

M 

PfailaU C (?A^)«l-71— ( 4 9) 
"-'max 

where Af max depends on q and N via Eq. (|4*2")l . Thus 
the necessary critical polarization is approximately de- 
termined via 

Pfail onc q ( £ c necc, M,N) = 1— ^— . (50) 

IVi- max 
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For N > 1. Eas. (|22|) and ([31]) imply that the left hand 
side is independent of N and q and these only enter into 
the problem via M max . Our strategy for determining the 
scaling behavior of e c nocc with respect to N is to con- 
sider fixed M/M max as N varies. This is equivalent to 
maintaining an ensemble size needed to sustain a fixed 
failure probability for the classical algorithm. In the typi- 
cal scenario M max m TT\f~N /4 and upon multiplying N by 
a factor 7, the ratio M/M max can only be held constant 
by multiplying M by a factor of yfiy. More explicitly, for 
any 7 > 0, Eq. (|48|) implies that 



Pfail < 



P {a «on c q (Ccn CCC (M,N),M,N) 



(51) 



and the approximation of Eq. ()37|) implies 

£ C necc(V7M, 7 7V) « 7 - 1 / 4 e c nocc (M, TV) 

or equivalcntly 

e c nccc(M, 7^) « 7~ 1/4 £c „ccc(M/V7, JV). (52) 

This relationship is evident in the numerical data; the 
plots of Fig. [1] illustrate an instance where 7 = 10 4 . 
The sufficient polarization satisfies 



Pfail < 



suir 



n 
1 



1 



M 



(53) 



Data generated via numerical solution of Eq. (|53p , using 
the incomplete beta function representation of Eq. (|36p , 
is illustrated in Fig. Q] A notable feature of the data is 
that for much of the region in which M max /2 ^ M ^ 
M maj(1 the necessary and sufficient critical polarizations 
are within an order of magnitude of each other; if all that 
is of interest is the typical order of magnitude of e c for a 
given N then it suffices to compute either e c nccc or e c su ff 
for M a moderately large fraction of M max - By this ar- 
gument one can attain an order of magnitude estimate of 
the critical polarization, which lies between the bound- 
ing critical polarizations. Additionally, for much of the 
range ^ M ^ M max , the numerical data indicates that 
the sufficient polarization is roughly constant. 

In the limiting case e c su g — > 1 as M — > M max . Also 
as M — > 1, p fail onc q (e c suff , M, N) -> (1 - e)/2 and thus 
Eq. (f53|) implies that 



2 

1 

n 



1 - 



4V./V 



1 



(54) 



in the typical scenario. The scaling of the sufficient criti- 
cal polarization is more complicated than that of the nec- 
essary critical polarization as a result of the additional 
factor of 1/n in Eq. (|5"3"|) . A similar line of reasoning to 
that for the necessary polarization gives 

Pfail onc^ec suff (y/^M, 7#). \W 7^0 



riog 2 Jvi 

[log 2 (lN)] 



Pfail < 



\e c suff (M,N),M,N) (55) 



Here the terms involving logarithms prohibit immediate 
use of Eq. (|37|) as was done for the necessary polariza- 
tions. However, when [log 2 7] <C [log 2 N] , 



suff (vW 7^0 



/ °c su: 



ff (M,JV), (56) 



since the factor involving logarithms on the right of 
Eq. (|55|) is approximately unity. Again this evident for 
the case of 7 = 10 4 as illustrated in Fig. [TJ 

An additional illustration of the scaling properties of 
the necessary and sufficient critical polarizations can be 
attained by computing these at values of M for which 
the classical algorithm returns a fixed success probabil- 
ity rate as varies. In order to succeed with probability, 
^ p ^ 1, the classical algorithm requires Q = \pN] or- 
acle queries. The ensemble quantum algorithm uses com- 
parable resources for M — \Q/q\ ~ Ap/n^/N w pM max 
ensemble members. Necessary and sufficient critical po- 
larizations can be computed using N and the value of M 
computed from it as just described. This method of fix- 
ing M is the same as that resulting in Eqs. (|5"T|) and (|56|) . 
Dropping the first argument, which is computed from the 
second and p, Eq. (|5ip is equivalent to 



-1/4 



■XN ) 



(57) 



where N is any fixed database size that is sufficiently 
large for the approximations to be valid. A similar 
conclusion applies to the sufficient polarization. For 
p — 0.90, numerically generated data is illustrated, using 
logarithmic scaling, in Fig. [2l 



10 



15 



20 



25 



30 



logi N 




logi 



FIG. 2: Critical polarization for success probability of 0.90 
Open squares indicate sufficient critical poiarization and solid 
dots necessary critical polarization. The straight lines are 
obtained by least squares fits. Similar plots are obtained for 
other values of p, the primary difference being that the lines 
are shifted vertically. 

Least squares linear fits for these data give 



log 10 £ C nccc = 0.078 - 0.250 log 10 N 
logi £c suff = 0.387 - 0.247 log 10 N 



which are consistent with Eqs. (|5T|) : the discrepancy in 
the coefficient of logi N for e c su g probably arising from 
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the factor of 1 jn in the definition of the sufficient polar- 
ization. 

We now briefly consider an alternative method of com- 
paring the two types of algorithm. Suppose that the 
ensemble algorithm uses fewer resources than the prob- 
abilistic sequential search in the sense that (qM) a = Q 
where a > 1. The fact that the classical succeeds with 
certainty when Q N — 1 means that, in the typical sce- 
nario, the ensemble implementation can only outperform 
its classical competitor when a < 2. It is easy to show 
that for the necessary polarization, ^ M ^ M max , 
where now M max = \ L N^ l '° 1 ' /q\ an( i f° r the typical sce- 
nario this leaves M max ps 4A^ 1 ' q ~ 1 / 2 ) /tt. The scaling ar- 
gument can be repeated by considering multiplication of 
N by a factor 7 > 0. A constant classical failure proba- 
bility is then maintained by considering M multiplied by 
a factor of -y^ 1 /"- 1 / 2 ). Repeating the argument involving 
the single bit failure probability for the ensemble quan- 
tum algorithm ultimately yields 



[l/4-l/(2a)] 



Z (N ) (58) 



with a similar inference for the sufficient polarization. 

Current room-temperature, solution state NMR using 
pseudopure preparation schemes @, [H, [25[ attains typi- 
cal polarizations of e ~ 10 -5 . The numerical data plotted 
in Fig. [TJ and the scaling of Eq. (f5"Tj) and (|56"j) suggest that 
these standard realizations would only gain an advantage 
for databases of size N > 10 22 . In this respect, the re- 
quirements on the ensemble for an advantageous realiza- 
tion of the Grover search algorithm are less demanding 
than those for the Deutsch- Jozsa algorithm [2(| , where it 
was found that e ^ 0.866 was necessary in the best case, 
regardless of the problem size. It should be noted that 
the standard number of oracle queries for such database 
sizes is of the order of 10 11 , which forms a lower bound 
for the number of pulses in sequence of gates; sufficient 
coherent control over such a number of pulses is evi- 
dently beyond current NMR experimental capabilities. 
Nevertheless, compared to the situation of the ensemble 
Deutsch- Jozsa algorithm, our results are promising since 
an increase of one order of magnitude in polarization will 
result advantages for searches on a database whose size 
decreases by four orders of magnitude. In fact, for a dra- 
matic increase in p olarization such as that produced by 
newer schemes |26l . |27| producing a large initial polariza- 
tion of e ~ 0.9 would be sufficient for modest database 
sizes where N < 10 6 . Also, to date, all solution state 
NMR realizations of Grover's algorithm [ll|, [r| [28| have 
involved databases for which N ^ 8 and ensemble sizes 
of M ~ 10 20 . Thus M > M max in these cases and, ac- 
cording to our criteria, the ensemble quantum algorithm 
uses greater resources than a classical sequential search 
which will locate the marked item with certainty. 



V. CONCLUSION 

In conclusion, we have extended methods for statis- 
tical comparisons of single bit output ensemble quan- 
tum algorithms to classical probabilistic competitors to 
the Grover search algorithm. Our results indicate that 
the polarization required for the ensemble quantum algo- 
rithm to outperform known classical competitors scales 
as N^ 1 / 4 where N is the database size and the ensem- 
ble size is scaled so as to fix the failure probability of 
the classical algorithm. In absolute terms, the required 
polarization is modest compared to that required for an 
advantageous ensemble realization of the Deutsch- Jozsa 
algorithm. 
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APPENDIX A: SUCCESS PROBABILITY FOR 
TWO QUBITS 

We aim to compute the probability with which qubit 
1 returns the correct value k\ times and qubit 2 returns 
the correct value &2 times. Consider the joint measure- 
ment outcomes on each ensemble member; denote a typ- 
ical outcome by X2X1 where the subscript indicates the 
qubit number. Let Iq denote the number of times that 
00 occurs, l\ the number of times that 01 occurs, I2 the 
number of times that 10 occurs, etc .... The probability 
with which this occurs is 



(Iq, I1J2, h) 



1-eY (\ 



h '1-eV 2 /l + 3e N h 



4 J V 4 / V 4 
where the multinomial symbol is 



(Al) 



Given values for k\ and k%, various combinations of 
lo,h,h,h are possible; all are subject to the following 
constraints: 



^07 h,h,h ^ 
^0 + ^1+^2 + h = M 
h + h = ki 
h + h = k 2 - 



(A3a) 
(A3b) 
(A3c) 
(A3d) 



Thus 

h = h- h (A4a) 

h = k 2 - h (A4b) 

l = M - ki - k 2 + l 3 (A4c) 

and together with Eq. (|A3a[) these give 

h^ki+k 2 - M (A5a) 

fci > h (A5b) 

fc 2 Z 3 (A5c) 

Z 3 > (A5d) 

I 



APPENDIX B: UPPER BOUND ON FAILURE 
PROBABILITY 

The inequality 

Pfaii aii q (£, M, N) < n Pfan oae \s, M, N) (Bl) 

is clearly true for n = 1. We adopt the following notation 
for the n qubit case: Pr (s„, . . . , s 3 , / 2 , si) is the proba- 
bility that we succeed in identifying bit 1 correctly, fail 
on bit 2, succeed on bit 3 and so on. For n — 2, 

Pfail ali q (E: M > = Pr (/a, /l) + ?r (/a, Si) + Pr (s 2 , /i) 

= Pr(/ 2 ,/i)+Pr(/ a ,ai)+Pr(a 3 ,/i) 
+ Pr(/ 2) /i)-Pr(/ 2 ,/ 1 ). 
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Eqs. (|A5b[) and (|A5c|) are equivalent to 

is^min^i.Aa). (A6) 

while Eqs. (|A5a|) and (|A5d[) are equivalent to 

l 3 > max (0, fci + fc 2 — M). (A7) 

Any value of ^ such that max (0, k\ + fc 2 — M) ^ Z3 ^ 
min (fcx , fc 2 ) together with Eqs. (|A4p yield fci and fc 2 ■ Thus 



(A8) 



I 

However, for bit 2, 

Pfiui o„e q (^ = Pr (/a, A) + Pr (/2, *i) (B2) 

with a similar result for bit 2. Thus 

Pfi.il aii q (^ M, N) = 2p iail on *(e, M, N) Pr {h, fx) 
^2p failonc q (£,Af,iV) 

since Pr (/ 2 , /1) ^ 0. The general result follows by induc- 
tion on the number of qubits. For n qubits, 



min (k± ) 

Pr(fc!,fc 2 )= {M-ki-ki + h,hi-hM-h,h)\ 

; 3 =max (0,ki+k 2 — M) 



4 

min (fei ,^2) 

4M 



53 (M - fci - fc 2 + i, ki - I, k 2 - I, l)\ (1 - e) M "' (1 + 3e)' 



tmax (0,fci+fe 2 — M) 
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Pfai\ aii q ( e ' N) = Pr (/„, any outcome on first n — 1) + Pr (s„, fail in some way on first n — 1) 

= j? fail one q (e, M, iV) + Pr (s„, fail in some way on first n — 1) 

= p fail one q (e, M, iV) + Pr (s n , fail in some way on first n — 1) 
+ Pr (/„, fail in some way on first n — 1) — Pr (/„, fail in some way on first n — 1) 

^ Pfail one q ( £ ' "^0 + ^ >r ( S ™' m SOme wa Y 011 fi rS t H — 1) 

+ Pr (/ n , fail in some way on first n — 1) 
< Pfiui o„e q (^ M, iV) + (n — l) Pfail onc q (e, M, N) 



; ail one 



l (e,M,iV) 



(B3) 



which proves the result. 



J 



APPENDIX C: PERFECT POLARIZATION BOUND FOR SINGLE BIT FAILURE PROBABILITY 



The single bit failure probability is 

Pfail one 



M 



2 V k 

fc=M min v 
1 M 



2 t-^i \ h 

fc=M-M mi „ + l 



M\ / 1 - e cff \ k ( 1 + e cS x A 



(CI) 



For e = 1 and q = q st d, 



giving 



Thus 



(C2) 



1 - Eeff < 



2 " 2(iV-l) 
1 + £ c ff < 1 



N - 1 



EeflE < I- 



(C3) which implies that the first term in Eq (jCljl satisfies 
I 



M\ ( 1 - e cff \ / 1 + e off x 



2 ^ \k I \ 2 I \ 2 J ^ 2 ^ 2 k (N - 

k=M„ v ' v 7 x 7 v 



1 


M 

£ 




M 






X 


{- 




M 




^ mill 


1 / 


M 



1 



[2{N~ l)] M »in 



2(iV-l) VMmin + V' 



* mm 



) ' Ji-h ' 



" 2 ^/ mi J [2(iV - 2iV - 3 l ^ J 
I 

since Stirling's approximation gives [20| 

V" 1 x /r^\ / M \ 12 



k =o - \M min J V ttM 



< C5 > * -^zb («) 
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for M ^> 1 and thus the first term in Eq (|C1[) is bounded 
by 



M 



2 ^ I fc 



fc=M min 



1 / 2 2(iV- 1) 



2 V ttA/ 2JV - 3 V TV— 1 



2 

M/2 



Af-fc 



(C7) 



A similar argument applies to the second term in 
Eq (|Clj) . giving the same bound. Thus 



2 2(N-1) 



nM 2N-3 \N — 1 



M/2 



(C8) 



where 



r(x)r( y ) y 



p = (1 - e)/2, and 



(D2) 



APPENDIX D: GAUSSIAN APPROXIMATION 
TO THE QUANTUM FAILURE PROBABILITY 

The single bit failure probability can be expressed ex- 
actly in terms of the incomplete beta function 

Pm onc q (e> M, N) =i [l p (M min , M - M min + 1) 

+ I p (M-M min + l,M min )] 
(Dl), 



M + 1 



Substituting from Eq. (jE)2|) into Eq. (|D1[) gives 



(D3) 



Pf ail o„e q (^^^)= o 



r(M+ 1) 



^-!(l-t) 



2 r(M min )r(M- M min + l ; jo 



Redefining the variable of integration via £' := 1 — 2t gives 

1 r(M + 1) 



M-M„ 



+t 



M-M„ 



M mill -11 



dt. 



ftailone^,^) 



2 m+i r(M min )T(M - M min + 1) 



/ [(l-i') Mm '"" 1 (l + M 



+ (l-t') A '~ A/min (!+*') 



'\M min — 11 At , 



1 



r(M + 1) 



However, 



2^+1 r(M min )r(M - Af min + 1) 

r(M + 1) 











■]df- f [. 






Jo 





Pfailone q (0,M,/V) = 



and thus 



P f ailo„e q (^^,iV) = -- 



1 1 



2 A ^+i T(M min )T(M - M min + 1) J 
T(M + 1) 



dt' 



2 2 M +! r(M roin )r(M - M n 



(D4) 



(D5) 



(D6) 



- / [(1-1/) M »>>*- I (1 + 1?) M - M "*» 
1) ./o 

+ (1 - t ') M - M —{l + t') A ^»-i] dt'. (D7) 



The gamma function satisfies T(n + 1) = n\ for integral For M ^> 1, M m ; n ss M/2 and Stirling's approximation 



n, giving 



T(M + 1) 



Ml 



r(M min )r(M - M min + 1) (M - A/ min )! (M min - 1)! 

(D8) 



gives |2C 



r(M + i) 



r(M min )r(M - M min + 1) 



(D9) 
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Considering the expression of Eq. (|D7[) . each of the expo- 
nents within the integral are approximately M/2. Thus 

Pmo^M,N) ^l-^Ej\i-t'^dt'. (D10) 

With t" := y/Mt', Eq. (|DlO| implies 

1 i rsVM / .//2\ 

(Dll) 

and for M 3> 1 the integrand approximates an exponen- 
tial function, giving 

P fail o„c q (^M,iV)«i--^ / e-/ 2 dt. (D12) 

V27T JO 

APPENDIX E: NECESSARY CRITICAL 
POLARIZATION VS M 

We show that e c n ccc(M, A/") is a monotonically increas- 
ing function of M . In general e c nocc satisfies 

Pfcfl onc q (^c nccc(M, AT), 9) iV)M = 1 - ^. (El) 



For A/ odd and any e it can be shown [2(| 

that Pfailono q (£,Af,A0 = Pfaaoiie q (e,M + 1,7V). 
The right side of Eq. (|Elj) decreases as M in- 
creases. Thus p fai] onc q (e c necc (M + 1, AT), M + 1, N) < 
Pfaiionc q ( £ c„ccc(M,iV),Af,7V) for odd M and since 

Pfaii onc q ( e ' M, N) decreases monotonically as e increases, 
this implies that, for odd M, e cnccc (AI + 1,N) > 
£c nocc(A-f, N). It remains to consider £ C nccc(A/ + 2, AT) 
versus e c necc (M, N) for odd M. Consider the ratio of the 
success probabilities, 



r(e, M, TV) := Psucccss onc £ ' '"' ' (E2) 

Success allW^) 



where p success aU c (gM, AT) = qM/N is the probability 
with which the classical sequential search succeeds 
and p succe ™ q (e>Af>Ar) = 1 - p M ono q (e, M, A 7 '). 
We shall prove that Ar(e, M, AT) := r(e, M + 
2, AT) - r(s,M,N) < whenever M is odd. 
Then r(e cnecc (M,N),M,N) = 1 implies that 
r(e c „ccc(M,A0,M+2,A0 < r(e cnccc (M,N),M,N) = 1. 
Thus n ecc(Af + 2, AT) > £ c nocc (A/, AT). Now, for A/ 
odd [U, 



success one 



(e,M,N)= 



/-_ M+l 

M + l M + l 



M-k 



where p = (1 + e)/2. In these terms 

Ar(e, A/, AT) 



AT 


1 


q 


M + 2 

"a7 /p 





M + 3 M + 3 



Then 



/M + 3 M + 3 



£(^_+3) /%(M+l)/2 (1 _ t) (M + l)/2 di 



r(^±i + i)r(^ti + i) 



and thus 



Ar(e, A/, N) = — ^± ^ f /_L_ ,(m+i)/2 (1 _ 4) (m+i)/2 _ J_ t (M+i)/2 (1 _ 4) (M + i)/ 2 l d 

-T-1/ M+1 M2 /„ I M + l M J 

4 



at r(M + i) r p 



q [T(M±1)]2 J Q [M+l 



t(l-t)-±\tt M -w{\-t) {M - 1)l2 &t. 



(E3) 



(E4) 



(M + 2)(M + 1) T(M + 1) fP (A/+1)/2 _ (M+1)/2 

[(M + i)/2p r(^i)r(^i)7o ( j ( j 



(E6) 



For ^ t ^ 1 it is easily verified that ^ t (1 — t) ^ Eq. (|E6|) is negative (except at the single point t = 1/2) 
1/4. Thus the term in parenthesis within the integral of 
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and Ar(e, M, N) < 0. 
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